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The interaction of shear bands with crystalline nanoprecipitates in Cu-Zr-based metallic glasses is 
investigated by a combination of high-resolution TEM imaging and molecular-dynamics computer 
simulations. Our results reveal different interaction mechanisms: Shear bands can dissolve precipi¬ 
tates, can wrap around crystalline obstacles, or can be blocked depending on the size and density 
of the precipitates. If the crystalline phase has a low yield strength, we also observe slip transfer 
through the precipitate. Based on the computational results and experimental findings, a qualitative 
mechanism map is proposed that categorizes the various processes as a function of the critical stress 
for dislocation nucleation, precipitate size, and distance. 
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I. INTRODUCTION 

Metallic glasses (MGs) have advantageous mechanical 
properties, such as a high yield strength and a large elas¬ 
tic limit, but suffer from brittle failure, especially un¬ 
der tension, at temperatures significantly below the glass 
transition.Under compression, improved ductility is 
found for composites of MGs and crystalline secondary 
phases, namely, for Gu-Zr-based,^“^ Cu-Ti-based,^ and 
Zr-Ti-based MGs.^“^^ Under tension, a small ductility 
with 1% to 2% strain is observed for Gu-Zr compos¬ 
ites containing nanocrystals.With a higher volume 
fraction of the crystalline phase, not only compressive 
but also significant tensile ductility is reported for Zr-Ti- 
based, Ti-based,^^ and Cu-Zr-based MGs.^^’^^ 

For dendritic precipitates, there is a correlation be¬ 
tween the location of dendrites and the occurrence of 
shear-band patterns.The improved ductility is gen¬ 
erally ascribed to the increased number of shear bands 
and their limited length given by the constraints of the 
crystalline phase.Thus, a high volume fraction of duc¬ 
tile crystalline phase improves the ductility in compres¬ 
sion and microindentation tests, while a brittle secondary 
phase does not.^^ This is confirmed by Song et al, who 
suggest a crystalline volume fraction between 40% and 
80% in Cu-Zr-based MGs for obtaining good mechanical 
properties,which is consistent with the fact that ten¬ 
sile ductility is observed only in glasses with high volume 
fractions of ductile crystalline phases. 

While the enhancement of macroscopic ductility of 
MGs with high volume fractions of a ductile crystalline 


phase can be explained by simple composite models, the 
influence of nanoprecipitates on the mechanical prop¬ 
erties of MGs containing a much lower volume frac¬ 
tion of crystalline matter is still not clear. Similar to 
the case of dendritic precipitates, shear-band patterns 
were also observed in glasses containing small spheri¬ 
cal crystallites with sizes around 2nm.^ These nano¬ 
precipitates can grow during deformation in certain 
metallic glasses.An increased growth rate of 
nanocrystallites in shear bands is observed,which has 
been related to enhanced atomic mobility inside shear 
bands.Deformation-grown nanocrystallites are ob¬ 
served to contain twins,which occur only in 
larger crystallites, e.g., with a size greater than 20 nm 
in a Gu-Zr-A1 MG.^^ These deformation-grown precipi¬ 
tates are the possible reason for strain hardening during 
nanoindentation, ^ as well as increased plastic strain dur¬ 
ing compression.^ It is proposed that the participation 
of the crystallites in the plastic deformation is the rea¬ 
son for the enhanced ductility: Wu et al. demonstrate 
that reducing the stacking-fault energy of B2 CuZr by 
alloying leads to increased twinning and higher ductility 
under tension.Pauly et al. propose that a martensitic 
transformation from the B2 phase to the B19' phase with 
a subsequent volume change is responsible for toughen¬ 
ing in Gu-Zr-based metallic glasses.This interpreta¬ 
tion, however, is not generally accepted. Gorteen et al, in 
contrast, note that the volume change of the martensitic 
transformation is very small and cannot contribute signif¬ 
icantly to toughening.They instead suggest that pre¬ 
cipitates increase plasticity by favoring the nucleation of 
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new shear bands over the growth of critical shear bands. 
Indeed, recent simulation and experimental results pro¬ 
vide evidence for the fact that crystal-glass interfaces 
serve as nucleation sites for shear bands and are there¬ 
fore responsible for the simultaneous nucleation of mul¬ 
tiple shear bands.In tensile tests and corresponding 
molecular dynamics (MD) simulations of nanolaminates 
of copper nanocrystals separated by thin Cu-Zr glass lay¬ 
ers, the crystal-glass interface acts as a source or sink for 
dislocations. Shear transformation zones (STZs) are ac¬ 
tivated by interactions with dislocations. 

Computational studies on the interaction of crystalline 
precipitates with shear bands provide further insights 
into nanoscale mechanisms. Lund and Schuh conduct 
quasi-2D molecular-statics simulations of a binary Len- 
nard-Jones system with a nanocrystal inclusion.They 
identify three mechanisms of deformation, depending on 
the ratio of shear-band thickness to crystal size. For small 
crystals, the deformation is accommodated either in the 
interface (for example, by rotation) or by dissolution of 
the crystal. For wide shear bands and intermediate crys¬ 
tal sizes, dislocations in the crystal nucleate at the in¬ 
terface. Finally, for crystals larger than the shear band, 
they observe homogeneous dislocation nucleation due to 



FIG. 1. An undeformed Zrss.sCuai.eAgr.oAlr.e sample after 
annealing in a TEM bright-held image. An overview of the 
sample is shown in (a), with clearly visible crystalline pre¬ 
cipitates. The inset (b) shows a magnihed view, indicating 
the transition from globular to dendritic morphologies during 
precipitate growth. The different brightness of parts of one 
precipitate in (b) is most likely due to the different geometri¬ 
cal orientations of different parts of the same precipitate with 
respect to the incident electron beam. The strong black-white 
contrast change observed for some nanocrystals in (a) is most 
likely due to twinning. 


stress building in the nanocrystal center. However, it is 
somewhat unclear how the observed homogeneous dis¬ 
location nucleation depends on the artificially induced 
shear band and the resulting stress state in the system. 
Shi and Falk conduct molecular-dynamics simulations on 
a monoatomic amorphous model system with a high frac¬ 
tion of bee nanocrystallites.They find that deforma¬ 
tion is induced at the interfaces and that shear bands 
bend around crystallites away from a direction of maxi¬ 
mum resolved shear stress. They also observe blocking of 
shear bands by crystallites. Because of the high fraction 
of crystalline phase, the system more closely resembles 
a nanocrystalline structure. The observation of the ini¬ 
tiation of plastic deformation at interfaces still matches 
the simulations by Albe et al?^ and underlines the im¬ 
portance of the crystal-glass interface in these composite 
systems. 

While it has been shown that interfaces promote shear- 
band nucleation and that precipitates can act as obsta¬ 
cles or can deform together with the matrix, there is 
no comprehensive study that investigates the influence 
of the size and number density of the precipitates. Fur¬ 
thermore, some mechanisms governing the interaction be¬ 
tween a propagating shear band and a preexisting pre¬ 
cipitate have been observed but not investigated and dis¬ 
cussed in detail. Therefore the goal of this study is to 
investigate the interaction of a shear band with preexist¬ 
ing precipitates in Cu-Zr-based MGs. In the experimen¬ 
tal part of the study, we anneal Zr-Cu-Ag-Al melt-spun 
ribbons to induce the formation of nanocrystalline pre¬ 
cipitates. We present transmission electron microscopy 
(TEM) images of the samples before and after defor¬ 
mation by cold rolling and identify the effects of crys¬ 
talline precipitates on the shear-band propagation. Us¬ 
ing molecular-dynamics computer simulations, we model 
composite systems with a metallic-glass matrix and crys¬ 
talline precipitates. We control the initiation of a shear 
band using a stress concentrator and put precipitates in 
its propagation path. The focus of the simulations is on 
the size effects of “hard” precipitates that do not par¬ 
take in the plastic deformation. Additionally, we study 
the shear-band interaction with “soft,” plastically de¬ 
formable precipitates. Finally, we derive a deformation 
map from the combined observations of simulations and 
experiments that classifies the observed mechanisms. 


II. EXPERIMENT 
A. Experimental setup 

We prepare metallic-glass samples of nominal compo¬ 
sition Zr 53 . 8 Cu 31 . 6 Ag 7 . 0 Al 7.6 from pure components (Cu: 
99.999%, Zr: 99.998%, Ag: 99.999%, Al: 99.999%; all in 
at.%) by prealloying using arc melting. After repeated arc 
melting with intermittent turning of the specimen to en¬ 
hance homogenization, the entire ingots are inserted into 
quartz-glass crucibles for melt spinning. The weight loss 
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FIG. 2. TEM image of a shear band in a Zr-Cu-Ag-Al ribbon 
after deformation. The green arrows mark the positions of 
crystalline precipitates. 

during alloying is minimal and subsequent composition 
analyses by energy-dispersive x-ray diffraction (EDX) 
confirm that the composition of the material is equal 
to the nominal composition within the accuracy of the 
measurement. For melt spinning, the ingots are induc¬ 
tively melted under an Ar atmosphere and the melt is 
ejected onto a rotating Cu wheel (tangential wheel ve¬ 
locity: 30m/s), resulting in completely amorphous thin 
ribbon samples of approximately 80-qm thickness. X-ray 
diffraction on the as-quenched ribbon samples does not 
indicate the presence of any crystalline fraction exceed¬ 
ing the sensitivity threshold of this method. Parts of the 
ribbon samples are cut to perform differential scanning 
calorimetry (DSC) measurements in a Perkin Elmer Dia¬ 
mond DSC device. Both isochronal and isothermal mea¬ 
surements under a purified Ar gas flow are conducted and 
in conjunction with microstructure analyses the time and 
temperature dependence of the evolving crystalline frac¬ 
tion is determined. On the basis of these results, the sam¬ 
ples for deformation processing are annealed in the DSC 
device at 410 °C for 3h. This thermal treatment results 
in the formation of nanocrystalline precipitates with an 
average diameter of about 70 nm and a number density 
on the order of 10^^ m“^. The resulting microstructure 
is shown in Fig. 1 . Given that the width of shear bands 
is on the order of 10 nm,^^ there is no straight path for 
propagating shear bands to avoid the interaction with 
nanoprecipitates in these samples. 

After the DSC heat treatment, the partially crystalline 
material is deformed by cold rolling at room tempera¬ 
ture in one step to true strain values of about 5 = 5% 
at room temperature, applying a strain rate of the or¬ 
der of 5 = ls“^. Subsequently, specimens for TEM in¬ 
vestigations are prepared by grinding, dimpling, and fi¬ 
nally precision ion polishing (PIPS, Gatan), using a low 
acceleration voltage of 2.5 kV and low incidence angles 
(< 4°) to minimize damage by the preparation process. 
The electron transparent samples are then analyzed in 


dark-field, bright-held, and high-resolution transmission 
electron microscopy modes in a Zeiss Libra 200FE TEM 
operated at an acceleration voltage of 200 kV. Since in¬ 
trinsic shear bands formed upon cold rolling could not 
be found during the TEM inspections of the thin foil re¬ 
gions, the shear bands generated as a result of sample 
preparation and/or subsequent sample manipulation are 
studied instead. Thus, the hnal sample state resembles 
the state in an in situ TEM experiment comparable to 
the work in Ref. 23, where the shear bands are generated 
at crack tips in the thin TEM foil. 

B. Experimental results 

Figure 2 shows a shear band in a sample after TEM 
preparation. It is noticeable that the shear band switches 
propagation directions in the vicinity of precipitates. Be¬ 
tween precipitates 2 and 3, the shear band has an addi¬ 
tional bend. The propagation path, in general, suggests 
that the shear band is “attracted” to the precipitates, 
possibly due to a stress field resulting from the density 
change on crystallization, which explains the change of 
direction between precipitates 2 and 3. Because of the 
processing of the samples, it can be excluded that the re¬ 
peated bending of the shear band results from a change of 
the external stress state: The deformation by cold rolling 
is performed in a single step, and the observed shear 
bands were created during TEM preparation, resembling 
an in situ experiment. As the path change in the presence 
of precipitates is rather large (in Fig. 2 around 45°) and 
correlated to the position of the precipitates, it is most 
likely induced by the presence of the precipitates. The 
literature supports this, as shear bands in homogeneous 
metallic glasses (Cu-Zr-based or otherwise) are straight 
on the length scale presented here.^^“^^ It is therefore 
clear that the precipitates play a major role in influenc¬ 
ing shear-band propagation and thereby the macroscopic 
plastic deformation of the material. TEM images of a sec¬ 
ond shear band, shown in Fig. 3, shed more light on this 
interaction. The shear band interacts with two precipi¬ 
tates: A crack follows the path of the shear band up to 
the first precipitate and the shear band continues from 
the first to the second precipitate. Because of the vis¬ 
ible crack opening near the first precipitate, it can be 
ruled out that the shear band originates from the sec¬ 
ond precipitate. While the first precipitate is passed, the 
second precipitate, which is encountered centrally, stops 
the shear band. This can also be observed at the end of 
the shear band in Fig. 2 . A detailed analysis of the gray¬ 
scale intensity distribution of the high-resolution TEM 
(HRTEM) image in Fig. 3(e) indicates that the shear 
band changes its path slightly near the crystalline precip¬ 
itate but does not proceed further or shows slip transfer 
into the precipitate. Additionally, a smaller, shear-band¬ 
like region emerges almost perpendicular to the previous 
propagation direction (yellow arrow). This indicates ei¬ 
ther a shear-band deflection, or a nucleation of a new. 
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FIG. 3. TEM and HRTEM images of a shear band. In the up¬ 
per left corner, the electron transparent hole is visible, which 
stems from the sample preparation for TEM. (a-b) Bright- 
held images of the same area, using a different tilt angle. The 
shear band appears as a white stripe, with part of it already 
cracked (bright white), (c) Dark-field image of the same area: 
The nanoprecipitates are visible and marked by circles. The 
crack stops at precipitate I; the shear band continues from 
there to precipitate II. (d) The power spectrum of precipitate 
II along the (110) zone axis revealing super lattice refiections 
of the martensitic B19' structure, (e) HRTEM image of pre¬ 
cipitate II, showing the shear band stopping at the precipi¬ 
tate. The contrast in (c) and (e) is enhanced to improve the 
visibility of the precipitates and the shear band. 


perpendicular shear band. Based on the highly local na¬ 
ture of the intensity distribution and based on the com¬ 
parison of the contrast of other precipitates, preparation 
artifacts can be excluded. Thus, this type of interaction 
observed here is part of the intrinsic interaction mech¬ 
anism between precipitates and advancing shear bands. 
The power spectrum in Fig. 3(d) indicates a precipitate 
with B19' crystal structure,which is consistent with 
prior observations in Cu-Zr-based MGs.^^ 

It is clear that the interaction of a propagating shear 
band with a distribution of crystalline precipitates de¬ 
pends on the actual stress state near the shear-band 
tip, the already-accommodated stress by the shear-band 
propagation and the details of the local distribution 
of crystalline precipitates, their sizes, and the residual 
stresses in the glass matrix due to the formation of 
the precipitates. Depending on these factors, there are 
various explanations for the observed paths: The path 
changes are caused either by a deflection of the shear 
band, by the blocking and subsequent renucleation of a 
new shear band, or by several nascent shear bands grow¬ 
ing together. The unhindered “passing” of a precipitate 
can be explained either by a temporary path change of 
the shear band or by the participation of the precipi¬ 
tate in the plastic deformation. Because of the necessar¬ 
ily limited amount of data that can be obtained in the 
experiment, we undertake MD simulations to test the 
aforementioned hypotheses and investigate their relation 
to the sample geometry. 





FIG. 4. Schematic simulation setup. The picture shows a cut 
through the three-dimensional simulation box in the xz plane 
at 1/2. A notch is inserted to control the origin of the shear 
band (yellow). The spherical precipitate is shown in blue. The 
box has open boundaries in the x direction and is otherwise 
periodic. A constant strain rate is applied in the z direction. 
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FIG. 5. Snapshots of a simulation with a 30-nm CuZr precipitate. The shear band wraps around the precipitate and continues 
unhindered. The glass matrix is colored according to the atomic strain. The precipitate atoms are shown in blue if they appear 
in the B2 structure; no defects are visible. The left column shows a cut through the middle of the precipitate. On the right, all 
atoms with rji < 0.3 are deleted. A video version is provided in Video 1. 


III. SIMULATION 
A. Simulation setup and analysis 

In order to gain more insights into the nanoscale mech¬ 
anisms of shear-band interaction with precipitates, we 
perform a number of MD computer simulations, which 
allow for an “m situ"^ observation of shear-band propaga¬ 


tion. We use the software lammps^^ to quench metallic- 
glass samples, insert precipitates, and perform mechan¬ 
ical testing. The simulated metallic glass is a Cu 64 Zr 36 
alloy modeled with a Finnis-Sinclair-type potential by 
Mendelev et al.^^ Metallic-glass samples with dimensions 
10 X 10 X 10 nm^ and 20 x 20 x 20 nm^ are prepared by 
melting the material at 2000 K and subsequent quenching 
to 50 K with a cooling rate of 0.01 K/ps. 
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VIDEO 1. Simulation of a shear band wrapping around a 30- 
nm CuZr precipitate as shown in Fig. 5. 

We use the sample geometry illustrated in Fig. 4 to 
investigate the influence of preexisting precipitates on an 
approaching shear band. Size effects are studied by vary¬ 
ing the diameter d of the precipitates and by controlling 
the distance between periodic precipitate images by vary¬ 
ing the box width 1. Sample sizes are 120 nm x / x 60 nm, 
with I = 10 nm, 20 nm, 30 nm, and 40 nm. With a sin¬ 
gle precipitate per simulation box, this corresponds to 
number densities for the precipitates of 1.4 x 10^^/m^, 
6.9 X 10^^/m^, 4.6 x 10^^/m^, and 3.5 x 10^^/m^, re¬ 
spectively. The Cu-Zr glass samples are replicated to 
reach the desired box dimensions, and spherical precipi¬ 
tates with diameters from 3nm to 40 nm were inserted. 
For this, a hole is cut into the glass matrix with a size 
chosen to accommodate the precipitate without overlap¬ 
ping atoms. For the CuZr precipitate, we use the experi¬ 
mentally observed B2 structure.^^”"^^ The B19' structure, 
which was found in the precipitates in the experimental 
part of this paper, is a distortion of the B2 structure. 

A notch controls the origin of the shear band and makes 
sure that it always hits the precipitate. Periodic bound¬ 
ary conditions are applied in the y and z directions and 
open boundaries in the x direction. The resulting struc¬ 
ture is equilibrated at 50 K for 2 ns with a barostat at 
ambient pressure in periodic directions. After equilibra¬ 
tion, no long ranging stress field around the precipitate 
is left; any mismatches are accommodated by the glass 
during the interface creation. 

The resulting composite samples are deformed at 50 K 
under a constant tensile strain rate of 5 = 4 x 10^/s in the 
2 ; direction up to a total strain of at least 10%. The trajec¬ 
tories from equilibrated to fully deformed samples are an¬ 
alyzed to observe the shear-band propagation path. The 
shear band is identified using the von Mises local shear 
invariant as implemented in the visualization tool 
OviTO.^^ Atoms with a local shear greater than around 
0.2 are assigned to the shear band. To observe plastic 
deformation events in the crystalline phase, we perform 
atomic structure identification, which can identify crystal 


FIG. 6. Comparison of the Orowan mechanism (a) with a 
shear band wrapping around precipitates (b). The slip plane 
is shown in gray, the precipitates in blue, the dislocation as a 
black line, and the shear band as a red plane. The precipitates 
can stop a dislocation because it must remain in a defined 
slip plane, while the shear band can temporarily leave its slip 
plane and continue unhindered afterwards. 

structures, stacking faults, and other defects. 

B. Wrapping and blocking 

The first observed interaction mechanism between a 
propagating shear band and a preexisting precipitate is 
shown in Fig. 5 and Video 1. Here, the shear band wraps 
around the precipitate like a carpet moving over a small 
obstacle. The particle does not deform and simply moves 
along with one half of the glass matrix. We call this mech¬ 
anism the wrapping mechanism. This kind of athermal 
mechanism is virtually unknown in crystalline materi¬ 
als. The closest analog in a crystal—dislocation climb— 
is purely thermally activated. It is therefore instructive 
to compare the two material classes, as shown in Fig. 6. 
In a crystalline material, a dislocation moves on defined 
slip planes. A change of slip plane is possible only for 
the screw components of the dislocation,^^ is connected 
with a high energy barrier, and is usually observed only in 
stage III work hardening.Therefore, the Orowan mech¬ 
anism applies: The dislocation is bent around the obsta¬ 
cle and finally forms dislocation rings [Fig. 6(a)]. These 
rings may pile up, thereby hardening the material. In a 
metallic glass, as in any isotropic material, all slip direc¬ 
tions are equivalent. Only an applied external stress dif¬ 
ferentiates the directions. Under tensile stress, the planes 
oriented in 45° angles towards the tensile axis experience 
the highest resolved stress. An obstacle can be avoided 
simply by temporarily and locally changing the slip path, 
thereby wrapping around the obstacle [Fig. 6(b)]. De¬ 
pending on the precipitate distance, this wrapping mech¬ 
anism is observed for precipitates with diameters smaller 
than 25 nm to 35 nm in our simulations. 

An alternative mechanism appears for increasing di¬ 
ameters and decreasing distances between precipitates: 
The shear band is blocked by the precipitate. This causes 
the simultaneous nucleation of a second shear band per¬ 
pendicular to the first one on the opposite side of the pre- 






7 



FIG. 7. Snapshots of a simulation with a 37.5-nm CuZr-precipitate. The shear band is blocked by the precipitate, while a second 
shear band is immediately nucleated in another plane of high resolved shear stress. The glass matrix is colored according to 
the atomic strain. The precipitate atoms are shown in blue if they appear in the B2 structure; no defects are visible. The left 
column shows a cut through the middle of the precipitate. On the right, all atoms with 77 ^ < 0.3 are deleted. A video version is 
provided in Video 2. 


cipitate as shown in Fig. 7 and in Video 2. Again, there 
is no slip transfer to the crystalline phase. The paths of 
the shear bands in this simulation and in the HRTEM 
image of the experimental sample [Fig. 3(e)] are compa¬ 
rable. Both show what looks like a shear band that starts 
to wrap around the precipitate but does not propagate 


further. In the experiment, though, no fully formed shear 
band appears perpendicular to the original one. This may 
be a result of the more complex stress state or the fact 
that a new shear band can nucleate at another precipitate 
that is not visible in the images. Still, a region resembling 
a nascent shear band appears perpendicular to the orig- 





























VIDEO 2. Simulation of a shear band being blocked by a 
37.5 nm CuZr precipitate as shown in Fig. 7. 


inal propagation direction, strengthening the agreement 
with the simulation results. 

A systematic investigation of the parameters of pre¬ 
cipitate distance and size reveal a clear correlation with 
the mechanism. For a more quantitative analysis, we de¬ 
fine an empirical parameter A, which is given by the ra¬ 
tio of the cross-sectional area A of the precipitate di¬ 
vided by the distance of precipitate centers I as shown in 
Fig. 8(a). Similar to the derivation of the Orowan stress, 
the crystalline volume fraction / = i^precipitates/^ can 
be estimated from the average crystallite distance by the 
relation^ ^ 


7-^ 
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( 1 ) 


Using that, we can express A only in terms of volume 
fraction and precipitate geometry: 



dl2 ■ 


( 2 ) 


For nonoverlapping precipitates (d < /), it is A = ndP/A 
and A reduces to 


A = »jx/7. 


( 3 ) 


When the precipitates overlap {d > /), we reduce the 
area A to remove the overlapping circle segments. In the 
limit / ^ 0, A corresponds to A/I of an infinite cylin¬ 
der parallel to the shear-band front. Figure 8(b) shows 
a contour plot of the A parameter as a function of the 
precipitate diameter and distance. The data points in the 
plot represent results from our MD simulations, divided 
into those showing the wrapping and those showing the 
blocking mechanism. For the simulation geometry used 
in this work, there exists a given A that clearly separates 
the two mechanisms: 


Acrit ~ 12.65 nm 


( 4 ) 



FIG. 8. Influence of the precipitate size and distance on the 
mechanism, (a) Explanation of the parameter /, distance of 
the precipitates, and A, the area of a cut through the pre¬ 
cipitate. (b) shows a contour plot of the empirical geometry 
factor A = A/L The dashed line shows the critical value of 
A for a transition from the wrapping to the blocking mech¬ 
anism. Data points are MD simulations. The data points at 
I = Onm are infinite cylinders along the y axis, represent¬ 
ing the case of overlapping spheres with infinitesimally small 
distances between their centers. 


This Acrit is not universal, as, for example, the distance 
between precipitate and notch is not varied. A test sim¬ 
ulation finds that increasing the distance from the notch 
and therefore increasing shear-band length favors the 
wrapping mechanism. 

Using these formulas, we can also estimate A for the 
experimental results. Given a number density of precip¬ 
itates n = 10^^ m“^ and particle diameters of around 
70 nm, we obtain 

/ = = nU(35 nm)3 1,8% (5) 

V o 

Aexp = TT X 35nm x \/l.8% « 14.7nm > Acrit- (6) 

This is consistent with the fact that the shear bands 
that are observed in the glass samples are blocked by 
the precipitates. Still, the value is close to Acrit, which 
means that slightly smaller precipitates (or precipitates 
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FIG. 9. Stress-strain curves of samples which exhibit the 
wrapping or blocking mechanisms. The arrows indicate when 
the shear band hits the precipitate. 

that are not hit centrally, virtually decreasing the pref¬ 
actor d) may be susceptible to wrapping. 

Figure 9 shows examples of the stress-strain curves of 
samples that exhibit either the wrapping or the blocking 
mechanism. Neither a pronounced ductility nor signifi¬ 
cant strain hardening can be observed. This lack of strain 
hardening is in accordance with experimental data for 
tensile tests on Cu-Zr-based metallic glasses with crys¬ 
talline precipitates.It also fits with a recent study 
of compression tests of Cu-Zr-based metallic glass, that 
finds an effect of particle size on mechanical parameters 
but no particle hardening. 



FIG. 10. Mechanical dissolution of a 3-nm copper particle in 
a shear band. The gray-scale color coding shows atomic strain 
r]i from 0.0 (black) to 1.0 (white). The precipitate is shown in 
color: Light blue atoms are in the fee structure, green atoms 
are in a stacking fault, and dark blue atoms are disordered. 
The arrows indicate the shear direction. Subfigure (a) shows a 
cut through the middle of the nanoprecipitate, while (b) shows 
only the atoms that initially belonged to the nanoprecipitate. 



VIDEO 3. Simulation of a shear band cutting a 30-nm copper 
precipitate as shown in Fig. 12. 

C. Plastic deformation of the crystalline phase 

The precipitates discussed until now were all “hard,” 
i.e., not susceptible to plastic deformation under the 
simulation conditions: Because of the high antiphase¬ 
boundary energy in the B2 structure, superdislocations 
or twinning with respect to a martensitic transformation 
would be needed for a plastic deformation of the precip¬ 
itates. The stress available at the shear-band tip is not 
sufficient to nucleate these defects, which explains why 
no plastic deformation of the crystalline phase can be ob¬ 
served in our simulations. In experiments, deformation- 
grown precipitates show twinning defects which may 
cause softening effects and make the precipitate suscep¬ 
tible to plastic deformation.^^ As a model for a softer 
precipitate, we thus exchange the B2 crystal phase for 
fee copper. 

For precipitates that are small relative to the 
shear-band width, the nanocrystals undergo mechani¬ 
cal dissolution.This is also observed in our setup 
with 3-nm particles as shown in Fig. 10. In samples with 
larger diameters, the precipitates do not deform plasti¬ 
cally but instead show the same wrapping and blocking 
interactions as described before for the “hard” precipi¬ 
tates. Even if the (111) glide plane is oriented parallel 
to the shear-band direction to maximize resolved shear 
stress on the preferred fee slip plane, we do not observe 
slip transfer into the nanoprecipitate. The corresponding 
stress-strain curves are shown in Fig. 11(a). To explain 
this, we estimate the critical stress for heterogeneous dis¬ 
location nucleation in fee copper by shearing a nanowire 
on the (111) plane in the [110] direction [see Fig. 11(c) 
for the simulation setup]. For this, we hold the lower 
layer of atoms fixed and move the top layer of atoms 
with a constant velocity in the [110] direction to achieve 
volume-conserving shear. The diameter of the nanowire 
is 10nm, i.e., on the order of the smaller precipitates 
to maximize surface effects. The resulting shear stress 
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(a) engineering tensile strain e 

0 0.04 0.08 0.12 



(b) engineering tensile strain e 

0 0.04 0.08 0.12 



engineering shear strain 7 


(c) 



FIG. 11. Stress-strain curves for the composites with fee copper precipitates and for a 10-nm copper nanowire, (a) Curves 
from simulations with the Mendelev potential. The copper precipitates do not deform plastically and show similar behavior 
to the CuZr precipitates discussed earlier. The curve of the copper nanowire explains this: The yield stress, and therefore the 
critical shear stress, is higher than the highest resolved shear stress in the steady state in the composite. To be able to compare 
the tensile stress in the composite with the shear stress in the nanowire, the tensile stress axis is scaled with a factor 0.5 
corresponding to the Schmid factor for the plane of highest resolved shear stress, (b) Curves from simulations with the Ward 
potential. Now the critical stress for heterogeneous dislocation nucleation is lower than the steady-state stress in the composite 
and the precipitate deforms plastically, (c) The simulation setup for the nanowire. The wire is sheared in the [110] direction 
of the fee crystal structure on the (111) plane. The red atoms are fixed and the atoms on the top are shifted with a constant 
velocity to shear the nanowire. 


over shear curve is also plotted in Fig. 11(a). The yield 
stress of the nanowire, Tcrit, is an estimate for the stress 
needed for heterogeneous dislocation nucleation at the 
glass-precipitate interface. For comparing the shear stress 
in the nanowire with the tensile stress in the composite, 
we assume a Schmid factor of 0.5 and scale the tensile axis 
by a factor 0.5 compared to the shear stress axis. This 
graphically estimates an upper bound of the shear stress 
Tm that arrives at the precipitate, corresponding values 
are given in Table I. The upper bound of the estimated 
maximum resolved shear stress in the composite is com¬ 
parable to the lower bound for heterogeneous dislocation 
nucleation. While this suggests that dislocation nucle¬ 
ation may be possible, it is important to keep in mind 
that the nanowire shear test provides only a lower bound 
in an idealized case and that the actual shear stress avail¬ 
able to nucleate a dislocation may be lower than O.Scr^, 
which is why no plastic deformation of the precipitate 
is observed. Evidently, the value of Tcrit in the Mendelev 
potential is much too high, and therefore the Cu precip¬ 
itates are “harder” than expected, which is a deficiency 
of the potential model for pure Cu. 

Because of this, we switch to a different potential, 
which provides a better description of crystalline Cu. The 
Finnis-Sinclair-type potential by Ward et al.^^ is created 
by using preexisting potentials for the elemental phases^"^ 
and fitting the cross terms to the inter metallic phases. As 
shown in the Appendix, this potential has a more realis¬ 
tic unstable stacking-fault energy and critical stress for 
homogeneous dislocation nucleation than the Mendelev 


TABLE I. Comparison of maximum stress am and the cor¬ 
responding resolved shear stress Tm with the critical shear 
stress Tcrit for heterogeneous dislocation nucleation in copper 
in the Mendelev and Ward potentials. The calculation of the 
resolved shear stress assumes a Schmid factor of 0.5. 



am (GPa) 

Tm (GPa) 

"Tcrit (GPa) 

Mendelev 

3.2 

1.6 

1.55 

Ward 

2.9 

1.45 

0.75 


potential. As shown in Eig. 11(b) and Table I, the criti¬ 
cal stress for heterogeneous nucleation is much lower than 
even the steady-state stress in the composite, easily al¬ 
lowing plastic deformation of the particle. 

The results for a 30-nm copper precipitate are shown 
in Eig. 12 and in Video 3. The precipitate is cut by the 
shear band, and slip transfer through the particle can be 
observed. This mechanism replaces the previously dis¬ 
cussed blocking of the shear band if the nanoprecipitates 
are “soft”: Eor the crystal to partake in the plastic defor¬ 
mation, dislocation nucleation must be possible at shear 
stresses below the highest resolved shear stress in the 
metallic glass at yield. Despite the participation of the 
crystalline phase in the plastic deformation, the stress- 
strain curve in Eig. 11(b) shows the distinctive stress drop 
connected with a single critical shear band and no strain 
hardening. The reason is that in this setup the crystalline 
phase accounts only for roughly 5vol% of the sample. 
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FIG. 12. Snapshots of a simulation with a 30-nm copper precipitate. This simulation uses the Ward potential, in which the 
critical stress for dislocation nucleation is realistic. As the precipitate is sufficiently soft, the shear band can cut through it. 
The glass matrix is colored according to the atomic strain (the same scale as Figs. 5 and 7). Defects in the fee crystal structure 
are colored according to the legend. The left column shows a cut through the middle of the precipitate. On the right, all atoms 
with r]i < 0.3 and all fcc-coordinated atoms are deleted. A video version is provided in Video 3. 
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This means that the macroscopic mechanical properties 
are still dominated by the metallic glass. As the shear 
band can simply cut through the crystal, the precipitate 
poses no obstacle to the percolation of the critical shear 
band. For a larger crystalline volume, a ductile crystalline 
phase could possibly also constrain the shear bands. 


IV. DISCUSSION 

Using TEM imaging, we observe shear-band bending 
around or close to precipitates, an attraction of shear 
bands to the precipitates, and shear bands being blocked 
by precipitates. In our MD simulations, we find four 
mechanisms of interaction between shear bands and pre¬ 
cipitates: 

(i) precipitates that are small relative to the shear- 
band width dissolve mechanically, 

(ii) shear bands can wrap around precipitates, 

(iii) shear bands are blocked by precipitates, and 

(iv) shear bands cut through precipitates, and slip 
transfer into the crystalline phase takes place. 

Which of these mechanisms is active for a given precip¬ 
itate depends on the competition between the propaga¬ 
tion of the existing shear band, the heterogeneous nucle- 
ation of a new shear band, and the heterogeneous dis¬ 
location nucleation in the precipitate. The wrapping-to- 
blocking transition can be quantified by the parameter 
A = A/1 (X A^fj /d. Below A = Acnt, the wrapping mech¬ 
anism is favored. This value can be be explained by the 
following simple argument. When the shear band reaches 
a precipitate which does not deform, the shearing of the 
sample momentarily stops. The stress tsb in the shear 
band resulting from the externally applied tensile stress 
(Text amounts to 


(^) 

This stress acts mainly on the shear-band front, allowing 
us to write 


^ext ~ X I X hsB = I ^SB, (8) 

where hsB is the width of the shear band. At the moment 
that the shear band hits the precipitate, the force Text 
must be equal to a reaction force Fback from the precip¬ 
itate {actio est reactio). Using the projected precipitate 
area A (cf. Fig. 13), we can convert that force into a 
normal stress: 


a A — 


-Fback 

A/2 ■ 


( 9 ) 


We assume that Text predominantly acts on one half 
of the obstacle (area A/2), which is supported by the 



FIG. 13. Projection of the forces and stresses acting around 
the precipitate (blue circle) onto the xz plane. Shear bands 
are shown as hatched areas, where yellow signifies the arriving 
shear band, green the path for wrapping, and red the site for 
the nucleation of a new shear band. 


deformation pattern of the plastically deformed particle 
(Fig. 12 and Video 12). This back stress results in a shear 
stress Twrap ~ 0.5(7^ in the plane of the wrapping shear 
band (green shear band in Fig. 13). With Fback = ^ext, 
it is 


n -^ext ^ext ^ ^SB 
^ : zz — : — 


A/2 


A 


wrap • 


( 10 ) 


Twrap must surpass a critical value fwrap to allow the ini¬ 
tiation of the wrapping mechanism; otherwise the shear 
band simply stops propagating. This is not observed in 
our simulations, suggesting that (Text at yield is greater 
than 


<^ext 


2Af^ 


I hsB 


(11) 


The competing mechanism, blocking the shear band and 
nucleating a new one (red shear band in Fig. 13), can sim¬ 
ply be expressed by a critical shear stress fnud- Because 
of the low temperature in the simulation and a stress 
close to the yield stress, we consider only the athermal 
case and do not invoke a nucleation term which takes 
into account the relative volume of the interface. With 
^nuci = O.ScText: f^c transition from wrapping to blocking 
takes place where 

2fnucl = CText = ^ giving (12) 

I n-SB 

= (13) 

^ "^wrap 

This derivation also works in the case of externally ap¬ 
plied shear stress, by replacing (Text with 2rext- 

Assuming that nucleating a new shear band at the in¬ 
terface and propagating the wrapping shear band along 











13 


precipitate distance 



precipitate size 


FIG. 14. Schematic view of different mechanisms for the interaction of a shear band with a crystalline precipitate. With 
increasing precipitate sizes, the dissolution of the precipitate is first replaced by the wrapping mechanism. Depending on the 
critical stress for dislocation nucleation, wrapping is replaced by blocking or cutting. Wrapping can be favored by increasing 
the precipitate distance. 


the interface have similar critical stresses, we can simplify 
Eq. 13 to 

Acrit ~ ^SB- (14) 

Shear bands in Cu-Zr-based glasses have widths of 
around 10 nm,^^ which fits to the Acrit = 12.65 nm ob¬ 
served in our model systems. As stated earlier, the wrap¬ 
ping mechanism becomes more favorable again if the 
shear band is longer before it hits the precipitate. The 
reason is that this gives the shear band time to deviate 
slightly from its path, so that it does not hit the pre¬ 
cipitate centrally, thus effectively reducing A. In prac¬ 
tice, this is not a big problem, as the free shear-band 
length is constrained to approximately I anyway due to 
the distribution of precipitates in the sample. While this 
derivation is only approximate, it seems to be sufficient 
to explain the observed phenomena and guide future ef¬ 
forts in tuning the mechanical behavior of crystal-glass 
composites. Furthermore, it is easily possible to explain 
the fourth mechanism, a slip transfer into the crystalline 
phase. The critical stress for heterogeneous nucleation of 
a dislocation in the precipitate fdisi must be provided by 
the shear band via tsb = 0.5(Jext- If f^disi < f^nucb we can 
simply replace fnud by fdisi in Eqs. 12 and 13, thereby 
replacing the blocking mechanism with the plastic de¬ 
formation of the precipitate. A lowered fdisi also lowers 

-^crit • 

A simple deflection of the shear band is not observed 
and seems unlikely, as any deviation of the shear band 
from its path leads to a reduced resolved shear stress 
and thereby to a driving force to put it back “on track.” 


Contrary to the experiment, a change of shear-band path 
towards the precipitate was also not observed. Because 
the precipitates in our simulations are inserted artifi¬ 
cially and do not have a large stress field around them, 
this seems reasonable. In thermally grown precipitates, 
a stress field due to density mismatch between glass and 
crystal seems likely. Still, due to the geometry of the MD 
simulations, the shear band has two equivalent propa¬ 
gation pathways from the notch but always chooses the 
one leading towards the precipitate. This seems to be a 
weaker form of the attraction observed experimentally. 

With these results, we can attempt an explanation of 
the experimentally observed phenomena. First of all, the 
blocking of the shear band is a one-to-one correspondence 
between simulation and experiment. Comparing Fig. 3(e) 
with Fig. 7, we can see that the path of the shear band 
looks identical. The shear band wraps partly around the 
precipitate but is then blocked and does not propagate. 
Contrary to the simulation, no fully formed shear band 
but only a small shear-band-like region appears at the 
opposite crystal-glass interface. This may be a result 
either of the more complex stress state in the experi¬ 
ment or the fact that other precipitates are available at 
which the new shear band may nucleate. For the winding 
shear-band path, we can now exclude a simple deflec¬ 
tion as discussed above. A possible explanation would be 
the concurrent nucleation of nascent shear bands at the 
crystal-glass interfaces which grow together into a single 
mature shear band. While the interfaces are known to be 
nucleation sources for shear bands,our simulations 
show that the nucleation of a shear band at a stress con- 
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centrator like a notch or a crack always takes precedence 
to nucleation at interfaces or surfaces. The shear bands 
shown in the TEM images all originate from crack tips, 
making it unlikely that the shear band shown in Fig. 2 
consists of several concurrently nucleated shear bands. 
This leaves the explanation that this winding path is a 
series of subsequent blocking and renucleation events. 

The observed Acrit corresponds to precipitate diam¬ 
eters somewhere between 20 nm and 40 nm, depending 
on the interparticle distance. This critical diameter is on 
the order of magnitude reported in several experimen¬ 
tal studies for twinning in B2 crystals in Cu-Zr-Al-based 
metallic glasses of 20±5 nm.^^’^^ It also fits an experimen¬ 
tal work on Al-based glasses, where crystallites growing 
during deformation are sheared apart when they reach 
a critical size of about 10 nm.^^ Cu 5 oZr 45 Ti 5 metallic 
glasses exhibit a critical size of about 9nm for twinning 
of B2 precipitates.^^ The interparticle distances in these 
experiments are on the same order of magnitude as for 
our simulations. This supports an explanation of a tran¬ 
sition from wrapping to slip transfer in these systems. 

Figure 14 summarizes the competition between the dif¬ 
ferent mechanisms. Mechanical dissolution of the crys¬ 
talline particles occurs only if their size is comparable to 
the shear-band size.^^ With further increasing precipitate 
sizes, the shear band can still wrap around the obstacle 
until the size reaches a threshold value. This critical size 
also depends on the precipitate distance, as discussed be¬ 
fore, expressed in the parameter Acrit /). If wrapping is 
no longer possible, a precipitate which reacts only elas¬ 
tically to the applied stress will block the shear band. If 
the precipitate is susceptible to plastic deformation, slip 
transfer into the precipitate will take place. 

Concerning the mechanical performance of such in situ 
composites with crystalline precipitates that originate 
from nucleation and growth within the glass, the cur¬ 
rent results suggest that the discussed geometrical ef¬ 
fects serve to improve the macroscopic mechanical per¬ 
formance. None of the presented mechanisms seem to in¬ 
hibit the percolation of a critical shear band, yet, catas¬ 
trophic slip along a shear band leading to complete fail¬ 
ure is delayed in the case of winding shear bands due 
to the increase of the shear-band path length as well as 
the raised activation barriers for slip along shear bands 
that have a more complex topology. The wrapping mech¬ 
anism does not pose an obstacle to shear-band propaga¬ 
tion, but can be avoided by appropriate adjustment of 
the crystalline volume fraction and precipitate diameter. 
“Soft” precipitates additionally open possibilities to ad¬ 
just the plastic deformation by participating in it. Con¬ 
sistently, by increasing the volume fraction of the ductile 
crystalline phase, the constraints on shear-band propaga¬ 
tion can be increased, immediate failure can be prevented 
(cf. Ref. 11), and the composite displays macroscopic me¬ 
chanical behavior according to a mixing rule, further al¬ 
lowing one to tailor the properties. 


V. CONCLUSIONS 

Using TEM imaging, we observe shear-band bending 
around or close to precipitates, an attraction of shear 
bands to the precipitates, and shear bands being blocked 
by precipitates. MD simulations reveal that the shear- 
band bending is most likely the result of the subsequent 
blocking and renucleation of shear bands. Moreover, we 
identify shear bands wrapping around precipitates and 
slip transfer into the crystalline phase. By describing the 
competition between the critical stress for wrapping, the 
nucleation of a new shear band, and the nucleation of 
dislocations in the crystal, we could derive a mechanism 
map for metallic glasses with nanocrystalline precipi¬ 
tates. This detailed description of shear-band propaga¬ 
tion not only helps to understand the mechanical failure 
of these composites but also aids in tuning them. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 


displacement (A) 

FIG. 15. Generalized stacking-fault energies (top) and the re¬ 
sulting shear stresses (bottom) for fee copper. The DFT values 
for the stress curve are from Ogata et and the corre¬ 

sponding stacking-fault energies are approximated by using a 
numerical integration of the stress data. 
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Appendix: Generalized stacking-fault energy in fee 
copper 

The generalized stacking-fault energy in fee copper 
is calculated by using both Mendelev^^ and Ward^^’^"^ 
potentials. Density-functional theory (DET) calculations 
from Ogata et al.^^ using the same method are used for 
comparison. This is plotted in Eig. 15. The Ward poten¬ 
tial much more accurately describes the stacking-fault 
energy than the Mendelev potential, which has a criti¬ 
cal stress for homogeneous dislocation nucleation which 
is more than two times too high. 
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